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Abstract 

The binding of clusters of metal nanoparticles is partly electrostatic. We address difficulties 
in calculating the electrostatic energy when high charging energies limit the total charge to a 
single quantum, entailing unequal potentials on the particles. We show that the energy at small 
separation h has a strong logarithmic dependence on h. We give a general law for the strength of 
this logarithmic correction in terms of a) the energy at contact ignoring the charge quantization 
effects and b) an adjacency matrix specifying which spheres of the cluster are in contact and which 
is charged. We verify the theory by comparing the predicted energies for a tetrahedral cluster with 
an explicit numerical calculation. 



PACS numbers: 



I. INTRODUCTION 




FIG. 1: a: Transmission electron micrograph of experimental superlattice structure containing lead 
sulfate and dark colored palladium nanoparticles showing formation of regular palladium clusters 
as sketched in the colored inset. Scale bar is 20 nm. Reprinted by permission from Macmillan 
Publishers Ltd: from Ref. 3 , Figure lj, courtesy D. V. Talapin. 

b: Transmission electron micrograph of a dodecagonal quasicrystal superlattice self- assembled from 
Fe2 O3 nanocrystals and clustered dark-colored 5-nm gold nanocrystals. Reprinted by permission 
from Macmillan Publishers Ltd from RefP Figure 2b, courtesy D. V. Talapin. 

In self-assembled lattices of nanoparticles one often encounters clusters of metal particles^ 
as shown in Figure [I] The remarkable stability of these clusters was argued to depend partly 
on states of nonzero electric charge^. For particles of nanometer scale, such states are dom- 
inated by the quantization of charge. The energy to add a single electron to a particle 
becomes large on the scale of the thermal energy ksT, so that net charge on a particle is 
atypical. Thus any net charge on a cluster is necessarily unevenly distributed over its par- 
ticles. Still, a net charge on one particle must polarize the surrounding particles, producing 
electrostatic attraction. This contrasts with the macroscopic case in which the available 
charge would be shared amongst the particles, producing repulsion. It is of great interest to 
understand what types of clusters are favored under this simple and novel binding mecha- 
nism. Mutual electrostatic interactions between spherical conductors and with surfaces are 
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of interest in space environments^ and in scanning probe microscopj^. Merrill et aP explored 
the interactions among charged colloidal particles in clusters in solution. 

Unlike most interactions of small particles, this electrostatic interaction cannot be re- 
duced to a pairwise potential energy. Charge on one sphere induces polarization on each 
nearby sphere. This polarization induces further polarization in other spheres. Since their 
separation is not large compared to their radius, the polarization cannot be accurately de- 
scribed by a dipole approximation. Instead, all the spheres carry a polarization charge 
distribution that must be found self-consistently to minimize the electrostatic energy. It 
is not known what types of clusters would be favored by this novel multi-body interaction 
mechanism. Recently Moore^ has provided a multipole formalism for calculating this energy 
and has explored the energies of simple clusters . 

These polarization effects are nontrivial even for the case of two isolated spheres. Nu- 
merical solutions by A. RusselP, by Pisler et aP and by Kalinin et aP have been developed. 
A case of special interest is that of identical spheres of radius R bearing equal and opposite 
charge q at separation h. At small separation h <C R the charge becomes concentrated 
arbitrarily strongly near the contact point. This concentrated contact charge creates a log- 
arithmically singular mutual capacitance^ c(h) of the form c(h) — > \R\og(aR/h), where a 
is a numerical constant. The resulting electrostatic energy ^q 2 /c(h) reflects this singular be- 
havior. This divergent contact charge complicates the treatment of clusters of spheres with 
different charges. Moore's recent work on such clusters^ shows a non-regular dependence of 
the energy on separation. 

Here we investigate the implications of the singular contact charge for the electrostatic 
energy £{h) of clusters of conducting spheres % at small separation h when the total charge 
Q resides on only one sphere, as sketched in Figure [2j We contrast this energy to the 
simpler equipotential case where the charge Q is allowed to pass freely between the spheres. 
Here there is no contact charge, and the electrostatic energy S e (h) varies smoothly with h. 
However in the case of interest where only one sphere is charged, new behavior arises owing 
to the appearance of contact charge. It is necessary to characterize this new behavior in 
order to find the desired electrostatic energy when the separations h are small. We find 
that the energy at contact remains finite and equal to £ e (0), but it acquires a logarithmic 
correction in h: 

£(h)^£ e (0) [l+A/c(h) + ...]. (1) 
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FIG. 2: Sketch of a cluster whose electrostatic energy is to be calculated. Sphere 1 has a net charge 
Q; spheres 2, 3, and 4 have no net charge. Spheres 1, 2 and 3 are in near contact. Sphere 4 is in 
near contact with spheres 1 and 3 only. Near contacts have separation h much smaller than the 
sphere radii. Contact charges q between spheres 1 and 2, and spheres 1 and 4 are shown. Similar 
contact charges between spheres 1 and 3 and spheres 2 and 3 are hidden from view. Right. Same 
cluster with the charge Q free to migrate between spheres. Spheres are at the same potential and 
there is no contact charge. The regular tetrahedron treated in Figure [3] is obtained by moving 
sphere 4 so that it contacts all the other spheres. 

Remarkably, A depends only on the equipotential charge distribution and an adjacency 
matrix of the cluster considered. 

We begin by reviewing the origin of the singular c(h). Next we define a capacitance 
matrix C(h) that gives the proportionality between the charges Qi and the potentials V{. 
We then isolate the singular contribution to C(h) and thereby derive the result of Eq. [Tj 
To gauge the importance of the logarithmic correction in practice, we work out the case of 
a tetrahedral cluster. Finally we comment on experimental implications and tests. 

II. MUTUAL CAPACITANCE OF TWO SPHERES NEAR CONTACT 

For completeness we recall the origin of the logarithmic divergence of the mutual capac- 
itance of two neighboring spheres of radius R, bearing equal and opposite charges Q. The 
potential difference between the spheres is denoted V. In the limit h/R -C 1, the capaci- 
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tance is dominated by the adjacent sections of the two spheres. Since the curvature there is 
very small on the scale of h, we may find the capacitance from this region via the Derjaguin 
approximation^. This approximation treats the system as a set of concentric annular ring 
capacitors, neglecting the slopes of the surfaces within each ring. At lateral distance x from 
the central axis, the electric field E is evidently V/y(x), where y(x) is the gap thickness 
at position x. Thus the surface charge density cr(x) = E/(Atc) = V / (Airy(x)) . To find the 
charge Q, we integrate a: 

Q = J 2n x dxa(x) (2) 

We note that the local height y(x) is given by x 2 + (R — (y — h)/2) 2 = R 2 so that 2x dx + 
2(R ~{y- h)/2) {-\)dy = or for R > y, 

2x dx — > R dy, (3) 

and 

Q -> f aR 2nV Rdy/(47iy), (4) 

J h 

where a R is some upper cutoff of thickness where the Derjaguin approximation breaks 
down. Thus, 

Q^\v R f h R dy/y = \v R log(a R/h) (5) 

as claimed. The capacitance c{h) = Q/(2V) thus goes logarithmically to infinity as h — > 0. 
Numerically, one can use image charge methods^ to determine the explicit form of c(h) for 
small h: c(h) ~ \R\og(l.2QR/h). 

From this capacitance we can infer the energy needed to separate the contacting spheres 
with charges ±Q. At contact, the energy 8(0) is given by \Q 2 jC '. Since C — > oo, we have 
vanishingly small £ at contact. At infinite separation we have the full Coulomb self energy 
2 \Q 2 jR. Thus with equal and opposite charges the polarization of the spheres cancels 
virtually all the electrostatic energy of the separated spheres. 



III. EQUIPOTENTIAL CLUSTER 

We now consider a cluster of n spheres labeled by i. We denote the set of charges Qi 
by the vector Q. There is in general a linear relationship between the charges Q and the 
potentials on the spheres V of the form 

Q = CV, (6) 
5 



where C is an n x n symmetric matrix. To set the stage for the problem of interest, we first 
consider the simpler problem of an equipotential cluster at potential V e . It is convenient to 
define a "uniform vector" u = (1, 1, 1, ....1). Then in this equipotential case, the potentials 
have the form V = V e u. Likewise, the total charge Q is given by Q ■ u. The equipotential 
cluster has none of the contact charge considered above, and thus none of the non-regular 
dependence on separation h. That is, C{h)u may be replaced by C(0) -u. Thus the charges 
Qi also depend smoothly on h. For a given total charge Q, the potential V e is then given using 
Q = u ■ Q = V e u ■ Cu. Evidently the equipotential capacitance C e is simply Q/V e = u- Cu. 

IV. ASYMMETRICALLY CHARGED CLUSTER 

We now fix the charge Qi on each sphere. For definiteness we may consider the case where 
sphere 1 has charge Q and the others are uncharged. We seek the corresponding potentials 
V that create this charge distribution. From this case, one may infer the case of general 
Qi by superposition. In general the potentials Vi are now unequal. The linear relationship 
Q = CV now includes contact charge and hence a non-regular dependence on separation h. 

We note that two isolated spheres i and j at separation h much less than the radius have 
charges given by c(h)(Vi — Vj). For our cluster of spheres, we define the contact charge 
between two spheres as this qij. Thus the total contact charge qi on a sphere can be written 
J2j Qij = c(h)DV . Here the adjacency matrix Dy is -1 for contacting spheres i and j, while 
Djj is the number of spheres contacting sphere i. We note that for n connected spheres this 
D has at least one null vector, since any uniform potential (proportional to the uniform 
vector u) produces no contact charge. Likewise the range space of possible contact charges 
q has at least one constraint: the sum of all the g« i.e.,u ■ q\s just the sum of the q^, which 
vanish pairwise. Thus u ■ q = 0. 

Within these restrictions the contact matrix cD is invertible- s . That is, every q with 
u ■ q = has a corresponding potential v such that q = cDv and v • u — 0. The Appendix 
establishes this invertibility property. 

Not all the charge in the sphere cluster is contact charge. Thus in general for a given set 
of potentials V the total charge Q = CV has the form q + Q. That is 

CV = c{h)DV + Q (7) 
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— * — * 

Evidently Q can be written in the form CV; here C is the capacitance matrix for the non- 
contact charge. We suppose that only the contact charge q diverges for given V as h — > 0. 
That means C(h) is regular as h — >■ 0. Thus for h <C R we may replace C(h) by C(0). 

We now note that the singular behavior of contact charge implies that V becomes uniform 
as c(h) — > oo, independent of Q. 

Q = CV + c(/i)DV (8) 

When we let c — > oo for fixed Q, any nonuniform V has a nonvanishing DV, and thus 
c(h)~DV diverges. Such a divergence is incompatible with the fixed limiting values of Q and 
CV. Thus V cannot have a nonuniform part; that is V = V e u for some V e . Recalling the 
discussion of equipotential clusters above, V e is the potential of the cluster with the same 
total charge as the desired cluster, with the spheres all at the same potential. 



V. CORRECTION FOR NONINFINITE c{h) 

For large c(h) the departure of V from the uniform V e u must be small. Accordingly we 
express V = V e u + v, where v is a small correction. Using this form 

Q = C(V e u + v) + c(h)~Dv (9) 

In this expression Cu may be simplified. This quantity is simply the non-contact charge 
under a uniform potential V e . However in this equipotential situation all the charge is 
noncontact charge, so that Cu is simply Cu. Using this simplification, we obtain an implicit 
expression for v: 

X = Q - V e Cu = (Q + c(h)T>)v (10) 
We note that the charge vector X on the left has a vanishing total charge u • X: 

u-X = u-(Q-V e Cu) = Q-V e u-Cu (11) 

Since V e is defined by Q = V e u ■ Cu, the total charge in X vanishes as claimed. 



Eq. 10 expresses the desired v in terms of the regular matrices C and D, which have no 
singular dependence on h. In order to find v from this implicit expression, we must invert 
the matrix C + c(/i)D. We discuss the invertibility of D in the Appendix. There we show 
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that D _1 X is well-defined for any X with u-X — 0. Moreover, the desired matrix C + c(/i)D 
is invertible apart from exceptional c(h) values, as the Appendix also explains. 

We now consider the leading correction for small h, i.e., large c(h). In view of the last 
term in Eq. [9j v must be of order l/c(h). To order (l/c(/i))°, the equation reads. 

Q - V e Cu = c(h)~Dv (12) 

Recalling that Cu = Cu, we now obtain the potential correction v: 

v = {l/c{h))B-\Q - V e Cu) (13) 

Since v is now determined, the state of the system for small separation h is determined. 

As anticipated, the potential is nearly that of a cluster with only the total charge Q 
constrained. To determine the logarithmic correction owing to non-infinite c(h) it suffices to 
know a) the adjacency matrix D and b) the charge distribution for the equipotential cluster, 
namely Cu. 

As h increases, further corrections become important. Further calculation via this scheme 
would require CV for nonuniform Vs. Such extensions would be cumbersome to implement. 
However, even without carrying out such an expansion, the adjacency matrix D can be used 
to produce a smoother behavior in h. If one has a scheme for computing C(h) for some 
range of h, then one may remove the singular behavior induced by c(h) by constructing 
C(h) = C(h) — c(h)T). This smooth C(h) can be used to find quantities of interest for 
values of h where C(h) itself would be poorly converged. 

From this v the electrostatic energy may readily be found, as we now show. 

VI. ELECTROSTATIC ENERGY 

Given this expression for the potential vector V, we may find the electrostatic energy £ 
for the cluster with a single charged sphere. This £ can be expressed as | QiVi = \Q • V ■ 
For the case of interest, only sphere 1 is charged. For convenience we define the vector 
1 = (1, 0, 0, 0, ...0). Now the energy £ = \Q V ■ 1, that is, 

£= l 2 Q{V e + v-l). (14) 

To leading order this energy is simply \Q 2 /{u ■ Cu) or ^Q 2 /C e {0), Here C e (0) is clearly 
larger than the capacitance C\ of an individual sphere. Thus £ is smaller than that of the 
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separated spheres, viz. \Q 2 jC\. Indeed, for a compact cluster of spheres, the capacitance 
is proportional to the radius, and hence to the 1/3 power of the number of spheres. For 
a cluster of many spheres, the binding energy of the cluster approaches the electrostatic 
energy of the charged sphere, as in the two-sphere case with zero net charge treated above. 
To evaluate the first correction in v, we may express v in terms of the uniform potential 



V e using the asymptotic formula of Eq. 10 



v- f-»- (l/c(h)) 1 ■ D _1 (Q T-V e Cu) 
Recalling that V e = Q/C e (0), this yields 



Cu 
u ■ Cu 



(15) 



(16) 



2 C e (0) V c{h) 

This energy has the form announced in Eq. [lj evidently the part (1 + A/c(h)) in Eq. [T] 
is the quantity in (...) above. 

This expression shows that one may isolate the singular part of the electrostatic energy 
for a cluster of conducting spheres close to contact, using nonsingular quantities which can 
be readily computed numerically. The energy at h = but without conductance between 
spheres is the same as for the equipotential case where conductance is allowed. This means 
that the imposed distribution of charge among the spheres has no effect on the energy at 
h = 0. 

Increasing h necessarily increases the energy. Thus the correction in l/c(h) is necessarily 



positive (though this is not obvious from Eq. 16). As seen from Figure [3] below, this increase 
can depend strongly on h. In general it depends on which sphere is charged. One expects 
that the increase with h is strongest for spheres that have many adjacent neighbors. 



VII. NUMERICAL ILLUSTRATION 



In order to demonstrate the behavior of the formula in practice, we show a specific 
example: a regular tetrahedron of contacting spheres, with a charge Q on one of them. 
Since each sphere contacts all three others, the adjacency matrix D is given by 

3 -1 -1 -1 
-1 3 -1 -1 
-1 -1 3 -1 
-1 -1 -1 3 



D 



(17) 



9 







h/R 



0.05 



0.10 



0.15 



<o -0.10 
-0.20' 



FIG. 3: Binding energy of a tetrahedron of unit spheres where one has unit charge. Units^ are 
such that an isolated sphere has £ = \. Points are from numerical energy minimization with charge 
represented as point charges, as described in the text. Left two points used 1999 points per sphere; 



remaining points used 499 points per sphere. Curve was obtained from Eq 16 The £ e (0) needed in 
this equation was obtained from the method of Moore^ and verified by the point charge method. 
The curve has no adjustable parameters. This £ e (0) is shown by the heavy mark on the vertical 
axis. For 5nm diameter spheres^ in vacuum with a single electron charge, this binding energy is 
5.0 kT at room temperature. 

The equipotential charges Cu are evidently all equal. The capacitance C e (0) of the 
equipotential cluster is greater than that of a single sphere. We report its numerical value 
below. Evidently the ratio Cm/ (it • Cm) = (1/4,1/4,1/4,1/4), and 1 - Cm/(m • Cm) = 
(3/4, —1/4, —1/4, —1/4). In order to compute 1 ■ D _1 ( 1 — Cm/(m • Cm) ), we must find the 
vector w which solves Du> = (3/4,-1/4,-1/4,-1/4) with w ■ u — 0. The solution is 



Finally, 



w = (3/16, -1/16, -1/16, -1/16). 



1 • D _1 (l - Cu/(u-Cu) ) = l-w = 3/16. 



(18) 



(19) 
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To test the range of validity of Eq. 16 , we numerically evaluated the energy £ for various 
separations h. We placed iV point charges of variable strength over the surface of each 
sphere and numerically adjusted these charges to minimize the mutual Coulomb energy of 
the charges under the constraint that the total charge on sphere 1 should be 1 while the 
others are zero. Separately, we calculated this energy with only the total charge on all 
spheres constrained to be 1. This gave an electrostatic energy at contact of 0.281. The 
binding energy of the charged sphere to the others is thus evidently 0.281 — |. We confirmed 
this value separately using the multipole method of Moore^ 

The dependence of S on h is shown in Figure 1. The prediction of Eq. [16] is accurate at 
the ten percent level out to separations h of roughly 0.1R in this case. 



VIII. DISCUSSION 



The preceding sections have explored a peculiar type of Coulomb interaction arising from 
the charging constraints encountered at the spatial scales of nanoparticles. Below we note the 
limitations of our work and suggest experimental situations where the interaction discussed 
here might nevertheless be relevant. 

In order to demonstrate the specific features our mechanism, we have considered the 
simplest example that shows the necessary features. First, a cluster of spheres like those 
considered here has charge polarization extending beyond the induced dipoles normally 
considered. Second, any excess charge on the cluster is dominated by single electron charges 
residing on one or another of the particles. Given these two features, one should observe 
the singular dependence on separation h found above. Our main aim has been been to 
show the form of this singularity and how its magnitude may be predicted from simple 
geometric information. The relative electrostatic energies of different clusters are important 
for determining their relative abundance and stability. For real experimental situations 
the relative abundance of actual cluster shapes doubtless depends strongly on several other 
factors as well. 

In real clusters, it is artificial to assume a single charge on a particular particle; a number 
of charge distributions likely have significant probability. If a single charge is present, it may 
reside on any particle of a cluster that doesn't require an extra energy much higher than 
ksT. Thus in practice one may need to consider an average over several charge positions in 
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order to determine the stability of a given cluster shape. 

In real materials a net charge on a cluster is only created in combination with a counter- 
charge elsewhere. In situations where the charging energy is large on the scale of ksT, these 
countercharges are normally located close to their opposites. In addition one expects ambient 
charges of both signs; these ambient charges screen the potential arising from any assumed 
net charge. Either the countercharge or the screening effects would greatly alter the relative 
energies of different clusters. Yet these effects of other charges need not be crucial. As for 
the countercharge effect, one may sequester the countercharge so that it remains far from 
the cluster in question. One example of such sequestration is in the layered semiconductors 
used to make two-dimensional electron gases^. Here the countercharge consists of strong 
electron donors or acceptors held many nanometers away from the free charges of interest. 
The same effect is achieved by making the counterions reside in large colloidal particles or 
micelles^. Their large size entails a low surface potential and hence weak interaction with 
the charges of interest. In these same colloidal systems, the high charging energy assures 
that ambient charge is minimal so that screening is not important. 

Naturally real clusters like those of Figure [T] experience other forms of interaction unre- 
lated to net charge on the cluster. The organic coronas^ used to to stabilize the particles 
exert interparticle forces. So do steric interactions with other neighboring nanoparticles. 
Dispersion forces and solvent-specific chemical interactions are also present. In order to 
make reliable predictions of cluster shapes, one would need to add these conventional inter- 
actions to the charge-induced interactions considered here. 

Experimental consequences of our clustering mechanism could potentially be found in 
the binary lattices like Figure 1 that motivated our study. If our mechanism is important, 
one expects a) cluster shapes with lower electrostatic energy as calculated above should be 
relatively more prevalent, and b) particles with a thicker ligand layer should be less strongly 
bound but have greater preference for specific charge sites. Still, the number of competing 
effects that determine the specific cluster shapes precludes any decisive predictions. 

Other simpler systems give a brighter prospect for decisive predictions. One such system 
is a dilute dispersion of nanoparticles in a nonpolar solvent 12 . One may induce charge 
separation by adding large counterions to the dispersion^. Then any nanoparticles with 
a net charge will attract neutral nanoparticles via the mechanism described above. If the 
counterions are sufficiently large and distant, their effects can be made minor. Then one 
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expects to observe clusters with relative abundance dictated by the electrostatic binding 
energies described above. 



IX. CONCLUSION 



As shown above, the electrostatic energy of a cluster of spherical conductors has a novel 
form when one conductor is charged and their separations are small. In the limit of small 
separations the energy is finite, but the corrections to this limit are logarithmically singular. 
Thus for real clusters where the separation is nonzero, it is important to know the singular 
contribution. Both the limiting energy and the corrections can be expressed in terms of 
non-singular operations. It appears from our numerical example that these small separations 
can have a significant impact on the binding of the clusters. In certain situations as noted 
above, this distinctive form of binding could be significant in determining the prevalent 
cluster shapes. 

With these methods in hand, one may readily compare the electrostatic binding energies 
of a range of cluster shapes. Our work on these comparisons is in progress. 
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Appendix: Existence of inverses in Eq. 10 



Our general equation for the nonuniform part of the potential v requires that the matrix 
inversion implicit in Eq. [TO] is well-defined. As we have seen, the matrix D is not invertible 



as it stands. Thus the invertibility needed for Eq. 10 needs to be verified. We first discuss 



the invertibility of D for the restricted v's and Q's whose uniform part vanishes (u ■ v — 0) 



13 



and whose total charge vanishes (u ■ Q = 0). The D matrix is known in graph theory as 
the Laplacian Matrix^, and its invertibility properties are well established. Here we explain 
these properties in the present context. We then argue that C + c(h)D is also invertible for 



the purposes of Eq. 10 



To analyze D, we picture a circuit of consisting of n nodes connected by capacitors 
of capacitance c. Every contact point between two spheres i and j is represented by a 
capacitor connected between nodes % and j. Then the equations Q = cDV are simply the 
circuit conditions of equilibrium. The capacitance matrix of this network is cD. The desired 
invertibility property is that for any Q with u-Q = 0, there exists a v with u-v = such that 
Q = cD v. Circuit theory^ gives conditions such that the charges Q specify the potentials 
v. Below we verify this for any connected set of nodes. 

For simplicity we consider a basis set of Q where an arbitrary node k has charge 1 while 
node 1 has charge -1. We first establish the result for a chosen subset of the final set of 
contacts and then show that it remains true when new contacts are added. We choose the 
initial subset to be a "spanning tree"^ such that there is only one contact path between 
any two spheres. Then all spheres along the path between k and 1 have contact charges 
±1 and all the potential differences from sphere to sphere along the path k — > 1 are equal 
to — 1/c. The remaining contacts, not on the k — > 1 connecting path, have no potential 
difference. Thus any sphere % not on this connecting path has the same potential as its 
connected neighbors. Indeed, all the spheres % that connect to a sphere j on the k — > 1 
connecting path share its potential Vj. Thus for this case, there is a potential vector for 
any given charge vector. The above argument only defines potential differences, Thus it 
specifies V only up to an additive constant. To define a unique V, we specify further that 
the sum of the potentials must vanish. Following our notation in the main text, we denote 
this restricted set of potentials as v. This shows the desired result for the spanning tree 
network. 

From this base we may inductively construct the new potentials that obtain after a single 
contact is added, eg. between spheres i and j. Given a set of charges Q (summing to 
zero), we presume that these potentials x are determined in the given network without 
the i-j contact. We then show that the v are also determined after the new contact has 
been added. We denote the difference Xi — Xj = A. It is useful to define the potentials 
for a particular Q, namely, a charge Q' at i and —Q' at j. This charge necessarily resides 
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on the existing contacts of sphere i and sphere j. The corresponding set of potentials is 
proportional to Q'; we denote them as w Q'. We denote the potential difference between % 
and j, i.e.,(wi — wj)Q', as B Q' . 

We now add the new contact, and claim that the desired potential v is a linear combination 
of x and w. This potential must satisfy the network equation for the new contact as well 
as those of the pre-existing contacts. Both x and wQ' satisfy these equations for the pre- 
existing contacts. When the new contact is added, a new contact charge q appears on sphere 
i and — q appears on j. The desired v must satisfy q = c{vi — Vj). However, the new contact 
must not change the total charge on i or j. Thus the pre-existing contacts on node i must 
change their total charge by an amount Q' = —q. The new network with the added contact 
must have the same potentials as the pre-existing network with the added charges Q' . This 
Q' generates the potentials wQ' = —wq. Thus the net potential vector is x — wq. The 
potential difference between i and j is then A — Bq. The contact equation for the new 
contact then requires cq = A — Bq i.e., q = cA/(l + cB) . Moreover, the corresponding 
potential vector v = x — w q satisfies both the new and pre-existing contact equations with 
the original charges Q. This shows that the charges determine the potentials in the new 
network, as claimed. Since q in this system depends only on potential differences, we may 
readily add a constant to all the potentials to assure that u • v — 0. 

Our general expression for determining v depends on inverting the matrix C + c(/i)D. 
This expression is our decomposition of the full capacitance matrix C, which we presume 
to be invertible. In order to be useful, this inverse must be readily calculable for generic 
values of c(h). We now show that this is true provided C is invertible. We may re-arrange 
the equation (C + c(h)D) v = Q to form 

(l + (l/c)CD- 1 )Dt; = g (A.l) 

— * 

where I is the identity matrix. The matrix in (....) is invertible (for u • Q — 0) for generic 
values of c. To show this, we first note that CD -1 is invertible, since its two factors are. 
Further, since the two factors are symmetric matrices, the product CD -1 is also symmet- 
ric. Thus CD -1 has a simple eigenvector expansion with no zero eigenvalues. We denote 
these eigenvalues as Aj. The full expression (...) also has a full eigenvector expansion with 
eigenvalues 1 + \/c. Provided none of the A, are equal to — c, all these eigenvalues are also 
nonzero. Since 1+ (l/c)CD~ 1 has a complete set of nonzero eigenvalues, it too is invertible, 
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as claimed. Then we can determine v by 

v — D _1 (i + (i/ c )CD- 1 )" 1 g 



(A.2) 
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